Recurrence Time Statistics for Finite Size Intervals 
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We investigate the statistics of recurrences to finite size intervals for chaotic dynamical systems. 
We find that the typical distribution presents an exponential decay for almost all recurrence times 
except for a few short times affected by a kind of memory effect. We interpret this effect as 
being related to the unstable periodic orbits inside the interval. Although it is restricted to a 
few short times it changes the whole distribution of recurrences. We show that for systems with 
strong mixing properties the exponential decay converges to the Poissonian statistics when the 
width of the interval goes to zero. However, we alert that special attention to the size of the 
interval is required in order to guarantee that the short time memory effect is negligible when one 
is interested in numerically or experimentally calculated Poincare recurrence time statistics. 
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The recurrence of trajectories to a neighborhood 
of a region in the phase space can be used to an- 
alyze important properties of dynamical systems. 
When this tool is applied to experimental or nu- 
merical generated data the limit of infinitely small 
recurrence interval (the Poincare limit) is never 
achieved. In this article, we present a new ef- 
fect that appears due to the finite size of the re- 
currence interval and changes the exponential de- 
cay of the distribution of recurrence times. The 
results are analyzed for the logistic and Henon 
maps but are expected to apply to a large class 
of chaotic dynamical systems. 



I. INTRODUCTION 

Since it was settled, the Poincare recurrence theorem 
has been the source of a number of paradoxes relating 
reversible microscopic dynamics of a system on the one 
hand and irreversible macroscopic behavior of this sys- 
tem on the other hand. An answer to these paradoxes 
was given by Boltzmann, who adopted the law of big 
numbers (N — > oo, where N is the number of degrees of 
freedom of the system under study) and the recognition 
that the Poincare Recurrence Time (PRT) to a highly 
improbable initial condition is too large to be observed 
in times normally available. Boltzmann's point of view 
was recently restated by Lebowitz (Q) in contrast with a 
new point of view based on nonlinear dynamics. In this 
new approach, rather than the limit N — > oo, the cen- 
tral role is played by the sensitivity to initial conditions 
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added to the idea that the Poincare recurrence time do 
not need to be very large to lie beyond the observable 
range limit 

Besides its fundamental importance for classical statis- 
tical mechanics 30), PRT statistics has been used, in 
the recent years, as a tool for time series analysis in a va- 
riety of areas ranging from economics to plasma physics 
SHU;©, and as a way of stud ying trapping properties 
in Hamiltonian systems ficl Hit Tl2|) . which is an impor- 
tant feature for anomalous transport processes FLU) . 
The series of recurrence times itself has also been sub- 
ject of fractal analysis lfl5l ITrl ITW) . All these aspects 
have brought a renewed interest in the study of recur- 
rence time statistics. 

General results have shown that the exponential-one 
law, e , holds for the cumulative probability of first re- 
currence times (scaled by the mean first recurrence time) 
of transitive Markov chains l|18|) . hyperbolic dynamical 
systems such as axiom A diffeomorphisms 1H9I) and for 
systems verifying a strong mixing property 1)201) . Further 
improvement to the study of recurrences to finite size in- 
tervals has been given by Galves and Schmitt <j2lj) . They 
have computed an upper bound for the difference be- 
tween the cumulative probability of first recurrence times 
and the exponential-one law. Moreover, they have shown 
that the right scaling of the first recurrence times should 
include an extra factor besides the mean first recurrence 
time. This factor lies between two strictly positive con- 
stants independent of the recurrence interval. Recently 
these results have been extended to unimodal maps l|22l) . 

The Poissonian statistics, or its cumulative equivalent 
exponential-one law for scaled recurrence times, is de- 
duced in the limit of infinitely small recurrence interval. 
However, in many of the recent applications mentioned 
previously the recurrence times are obtained either from 
numerical simulations or experimental data. In these 
cases it is unavoidable to use a recurrence interval with 
a finite size. The size of the interval is chosen in order 
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to obtain a sufficient number of recurrences to build the 
recurrence time (RT) statistics. In this article we are 
interested in the RT statistics to finite size intervals of 
chaotic dynamical systems. 

We begin using simple basic concepts of combinato- 
rial analysis to deduce the statistics of recurrences to a 
given finite size interval for random processes and chaotic 
systems with strong mixing. This statistics applies not 
only for the first recurrence time (RT) but for all the r- 
th recurrence time. We obtain this statistics, which we 
call binomial-like distribution of RT, as a result of a sim- 
ple combinatorial analysis problem. This distribution is 
valid for every size of the return region. When the prob- 
ability of coming back to the return region is very small, 
the binomial-like RT statistics reduces to the Poissonian 
statistics, commonly observed for PRT problems in the 
literature (|| L23|) . Since we adopt a combinatorial 
approach to deduce this statistics, almost no information 
about the dynamics of the system, but its strong chaotic 
mixing property, can be obtained from it. Therefore, dy- 
namical properties show their signature when the recur- 
rence time statistics deviates from the former ones. One 
of these deviations is particularly important for Hamil- 
tonian systems and concerns with a power law tail for 
long recurrence times (|23f) . In this article, we study a 
type of deviation which is related with the presence of 
unstable periodic orbits inside of the recurrence interval. 
We call it short time memory effect. This deviation orig- 
inates the extra factor, which should multiply the mean 
recurrence time in order to give the right scaling of the 
series of recurrence times as considered by Galves and 
Schmitt l)2lf) . It appears when a finite recurrence inter- 
val is considered. Although this deviation is restricted 
to a few short recurrence times, it changes the whole 
statistics. Moreover, we would like to emphasize that 
our deviation is in agreement with the bounds estimated 
in the previous works. 

The article is structured as follows: in section [HI 
we present the deduction of the binomial-like statistics, 
which is exemplified by a Gaussian stochastic process. 
The statistics for chaotic dynamical systems (logistic and 
Henon map) are calculated numerically in section IIIII 
where the short time memory effect is clearly observed. 
In section IIVI we explore the origins of this effect and 
how it changes the RT statistics. Finally, in section IV! 
we summarize our conclusions of this article. 



II. BINOMIAL-LIKE DISTRIBUTION 

Let f : M — > M be a homeomorphism with an in- 
variant measure /z(M) = 1. Given a region I C M 
with n(I) > 0, the Poincare Recurrence Theorem asserts 
that a trajectory, having started inside /, returns to I 
infinitely many times. The time interval, Ti, between the 
i-th and the (i + r)-th return is what we refer to as the 
r-th recurrence time. This time interval is just one of 
an infinite sequence {Ti : i = 1,2, ... , oo}, and we are 



interested in the statistics of this sequence. 

For convenience, most of the calculations are made for 
unidimensional systems. In this case, the interval / is 
defined as I(X C ,S) = [X c — S, X c + 6], as illustrated in 
FIG. n with a Gaussian random time series. When we 
have small values of 5, and thus a small probability /J,(I), 
we are dealing with the Poincare recurrence time. 
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FIG. 1 First recurrence time to the interval I for a random 
time series with Gaussian distribution. 

This article concerns the discrete time case, where the 
system is observed at a constant sample rate r = 1. A few 
adjustments are needed for the continuous time case |3 

In order to obtain the statistics for the recurrence time, 
consider the following simple problem: Let ei and e 2 
be two mutually exclusive events. The event e\ occurs 
with the constant probability [i and e 2 with the constant 
probability (1 — fi). Consider now a sequence of Ti trials 
{Sk ■ k = 1, 2, . . . , Ti} where Sk — &\ or S% = e 2 . What is 
the probability of having r events of type e\ and (Ti — r) 
events of type e 2 with the constraint that the last trial 
results is an event of the type ei? This kind of problem 
is known in the literature of combinatorial analysis as 
Bernoulli trials ( "repeated independent trials for which 
there are only two possible outcomes with probabilities 
that remain the same throughout the trials" (26)). 

The answer for this problem is the following: the prob- 
ability of having r events e\ and (Ti — r) events e 2 is 
/i r (l — /i) Ti_r . The last event must be of the type e\, 
then there are 

{Ti 1)! 



(T,-r)!(r-l)! 



ways of having (r — 1) events e\ in the previous (Ti — 
1) trials. Combining these results and suppressing the 
index i, since T t is just one of an infinite sequence of r-th 
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recurrence times, we have: 



P(T;r, M ) 



(T-l)! 



(T-r)i(r- 1)! 



M r (l-M) 



T-r 



(1) 



For a dynamical system with an invariant ergodic mea- 
sure and for which each step is independent from the pre- 
vious ones, it is easy to see that Eq. Q gives the proba- 
bility of r-th recurrence time if we consider the following 
analogy: between the i-th and the (i + r)-th return to 
the interval I{X C , 5) the trajectory spends Ti steps, each 
step (one trial) has a probability ^ = fi(I(X c , 5)) of be- 
ing in the interval I{X C ,8) (event e\) and a probability 
[1 — n{I{X c , 8))] of being outside I(X C , S) (event e2). 

Usually one is interested in the first recurrence time 
statistics (r = 1). In this case, Eq. Q gives 



P(T;1 )M ) = M (1- M ) 



T-l 



(2) 



that can be rewritten as 



P(T;l, M ) = -^e 1 "( 1 -^, 
1 fi 

which, by its turn, reduces to the Poissonian statistics 
P(T;l, M ) = M e^ T , (3) 

when /i — > 0. 

This statistics is the one commonly encountered for 
Poincare recurrences in chaotic dynamical systems JfJ 
filil) . The small fi(I(X Cl 5)) condition is, usually, satisfied 
when we take small values of 5. 

For any r, in the usual limit lf26l) . fi — > and T> 1, 
the binomial-like distribution (Eq. QJ) reduces to the 
Poisson-like distribution 



P(T;r,M) = ^&^. 
(r-l)\ 



(4) 



Eqs. (j2J) and ((HJ), as well as any other distribution of 
first RT, must satisfy two conditions: the first one is, 
obviously, the normalization, 



£P(T; 1, //) = !, 



(5) 



and the second one is the Kac's lemma 



24: 27): 



(T) = ^TP(T;1, M ) 



T=l 



(6) 



Although Kac's lemma is usually applied to closed Hamil- 
tonian systems, its original derivation l)24|) was based in 
general assumptions that cover a large class of dynamical 
systems including those we shall discuss in this work. 

The above conditions will be used in section llVl to take 
into account dynamical effects on the recurrence time 
statistics. 



To finalize this section we shall exemplify the binomial- 
like distribution obtaining the RT statistics of a stochas- 
tic process whose variable, x, is governed by a Gaussian 
density of probability: 



Pg(x) 



y/2na 2 



2a 2 



(7) 



and a = 0.2 (see FIG-EJ 
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FIG. 2 Gaussian density of probability with (a;) = and 
a — 0.2. The probabilities of returning to the regions Ji and 
I2 is given by fii and {12, respectively. 

The probability, fi(I(X c , 6)), of returning to an interval 
I{X C , 5) is given by: 



fi(I(X c ,5)) 



X c +8 



p G (x)dx, 



(8) 



x c -s 



in particular, for the intervals l\ and I2, shown in FIG.|2 
we have {X cl = 0.3, <5i = 0.02) and {X c2 = -0.2, Si = 
0.06), respectively. Their corresponding probabilities are 
(j,(h) = 0.02596 and /i(7 2 ) = 0.09679. 

FIG. El shows the statistics of recurrences to the in- 
tervals I\ and I2 for the Gaussian random time series. 
The solid lines correspond to the analytical results given 
by Eq. <^Q) and they are in a good agreement with the 
numerical results. 



III. CHAOTIC DYNAMICAL SYSTEMS 

Let's see now what happens when we consider the 
statistics of first recurrence of two deterministic dynam- 
ical systems. 

The first system we analyze is the logistic map, 



X n J r \ bx n {\ Xn) 



(9) 
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FIG. 3 The r-th RT statistics for a Gaussian random walk. 
We chose arbitrary values of r and the two intervals illustrated 
in FIG. [21 The solid lines are given by Eq. Q. 



in the completely chaotic regime (parameter 6 = 4). Its 
invariant probability density, 



Pl(x) 



1 



7T-\/ X(l — x) 



(10) 



is shown in FIG.0] In this figure, we also see the return 
interval I\ (X cl = 0.9, Si — 0.01) with measure — 
0.02124. 
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FIG. 4 Invariant probability density for the logistic map. 
The marked area, fii, represents the measure of the interval 
h(X cl =0.9,5i =0.01). 

Using the logistic map, we create a trajectory with 
N = 10 s points and look for recurrences to the interval 
h = [X c i = 0.9, tfi = 0.01], indicated in FIG. H The 
statistics of first recurrence times to this interval is shown 
in FIG. El 

As our second example of dynamical system we choose 
the Henon map, 



x n+1 = l-ax r 
y n +i — bx n . 



Vn 



(11) 



FIG. 5 Recurrence time distribution for a logistic map 
(ii = [X c = 0.9,5 = 0.01] and N = 10 8 points) trajec- 
tory. The dashed line is the binomiaHike distribution for 
fi(Ii) — 0.02124 (which coincides with the corresponding 
Poissonian distribution). The full line is the curve calculated 
to satisfy © and @ after T = n* = 10. The inset is an 
amplification of the small return times region and the squares 
are the analytic calculation of the return time probability as 
explained in FIG. [7] 



with the parameters a — 1.4 and b = 0.3. 
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FIG. 6 Recurrence time statistics for the Henon map with 
J(X C ,(5) given by X c = (0.304,0.210) and 6 = 10" 2 . The 
total of 10 7 recurrence events were considered. Dashed line 
is the binomial-like distribution for = (T) -1 = 0.00785. 
The inset is an amplification of the short recurrence times 
region. 

The distribution of first recurrence to a finite size in- 
terval for the Henon map is shown in FIG. In order 
to obtain this result, we eliminate the transient by iter- 
ating the Henon map 10 3 times and take the final point, 
X c = (x c ,y c ), as the center of the phase space interval, 
J, with radius 5 = 10 -2 . An initial condition inside the 
interval / is chosen and the Henon map is iterated until 
we obtain 10 7 recurrence events. 
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Both maps present distributions of recurrence times 
which fall exponentially (a straight line for a linear-log 
graphic) after some T > n* (n* w 10 for the logistic map 
and n* » 40 for the Henon map). What is indicated 
through the solid lines in FIG. El and FIG. El for the lo- 
gistic map and Henon map, respectively. Nevertheless, 
these solid lines do not coincide with the ones (dashed 
lines) given by the binomial- like statistics, Eq. J3J), which 
is indistinguishable from what would be the correct Pois- 

sonian statistics with (T) = for the particular size 

of the intervals used. In particular, the difference be- 
tween the right mean recurrence time and the inverse of 
the slope of the solid line in figure El is approximately 
9.3%. This difference indicates that an extra factor, be- 
sides the mean recurrence time, should be considered in 
order to obtain the right scaling of the series of recurrence 
times Q. 

The reason for this deviation is the following: In the 
previous section, we showed that the complete indepen- 
dence of iteration and the existence of an invariant mea- 
sure were necessary to obtain the binomial-like statistics 
of recurrence to an interval of finite size. However, de- 
terministic dynamical systems hardly fulfill the condition 
that an iteration is completely independent from the pre- 
vious ones. Therefore, dynamical systems have a kind of 
memory which affects the distribution of short recurrence 
times as we shall show in the next section. 
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FIG. 7 Analytic calculation for the return time probability 
for T — 2 in the logistic map. The region R2 represents the 
points inside / that return to I after two iterations of the 
map. We also see that Ri = 0. 



IV. MEMORY EFFECTS ON THE RT DISTRIBUTION 

In this section, we show how the memory effect changes 
the probability of short recurrence times, and how the 
latests are responsible for the deviation of the whole dis- 
tribution. 

We take the advantage of working with a simple dy- 
namical system, the logistic map, which gives us the pos- 
sibility of calculating analytically the probability of short 
recurrence times. 



A. Short Recurrences and Lack of Memory 

The procedure to obtain the probability of short recur- 
rence times, illustrated in FIG.0 consists in identifying 
the region, R n C /, of initial points whose trajectories 
return to I after n iterations of the map f(x), what 
is equivalent to the first recurrence of the map f n (x). 
Since we are specially interested in the first return time, 
we note that, if R n n R m 7^ for m < n, the points 
from the intersection will return to I in m, not n, it- 
erations. To avoid this error, we define a new interval 
R n = R n — {R n n R m ) for every m < n. The probability 
of a first recurrence time T = n is, then, given by: 



P 1 (T = n) 



V(Rn) 



(12) 



This method provides the probability of any first re- 
currence time. Nevertheless, the determination of the 
regions R n becomes a cumbersome task as n grows. 

Although the above procedure is not useful for calcu- 
lating the probabilities of first return when n becomes 
large, it is useful for showing how the chaotic dynamics 
simulates a random process for large recurrence times. 

When n grows, the region R n becomes the union of a 
large number, N, of disjoint sub-regions in the interval /, 
that is: 



N 



R n = \jRi, with fi(R n )=J2»(K 



(13) 



The more sub-regions we have, the smaller they are. 

In particular, as the map f n (x) is a polynomial of order 
2 n there are 2™ unstable periodic orbits, each one with its 
associated sub-region. The 2™ sub-regions are distributed 
in the interval [0, 1] according to the invariant density, 
Pl(x), given by Eq. (jTU|) . The sub-regions R l n are just 
the ones which are inside the interval /. 

For large n, there are N « 2 n /j,(I) sub-regions in the 
interval /, each one with the same measure, 



(14) 



The probability of having a recurrence to I in the time 



6 



T 



n is: 



(15) 



Putting Eqs. (jT3|) and itT^l in Eq. ijT5|) . we obtain the 
recurrence probability, P{T = n) i=a /i(/), independent 
of ri for large n, which is the same hypothesis assumed 
for the Bernoulli trials problem in section This ar- 
gument justifies the exponential decay of the recurrence 
time statistics after the short times. It must be empha- 
sized, however, that P(T) given by Eq. (|15|l is not the 
probability of first recurrence to / in the time T = n. 
This one is given by Eq. ljT2T) . where R n , instead of i?„, 
is used. It is in the calculations of R n that the memory 
effect appears, since the approximation N 2 n /j,(I) is 
not valid for n < n* . The same arguments hold for a 
general hyperbolic chaotic system, since the number of 
recurrence sub-regions is equal to the number of peri- 
odic orbits, and the latest grows as e hT , where h is the 
topological entropy. 



B. The Fitting of the Memoryless Exponential 

Let's see how the deviation of the shortest recurrence 
times affects the whole distribution. As it was shown 
in FIG. for greater values of T (after the decay of 
the memory, T > n*) the recurrence time distribution 
approaches a straight line in a linear-log graphics, what 
can be generally represented by an exponential, 



Pex P (T) = goe 



- 7 T 



(16) 



Since we know the mean recurrence time (or the measure 
of the recurrence interval), the above two parameters, go 
and 7, can be analytically obtained by using the two con- 
ditions presented in section^] namely, the normalization: 



n*-l 



T=l T=n* 

and the Kac's lemma: 

n* — 1 oo 



(T)= ]T TP 1 (T)+ TP exp {T) 



(17) 



(18) 



where Pi(T) is given by Eq. (|12|l or obtained directly 
from the data. The agreement of the exponential 
(Eq. obtained by this procedure with the linear part 
(in the linear-log representation) of the RT distribution 
for the logistic map, as shown by the solid lines in FIG.[S] 
and FIG. |SJ was verified for different recurrence intervals. 



C. Dependence of the memory effect on the interval 

In this sub-section, we explore the dependence of the 
memory effect on the position X c and size of the interval 



with a special interest in small intervals. As we argued 
before, only a few number of short recurrence times di- 
verge from a straight line in the linear-log representation, 
but their effects on the whole distribution is consider- 
able. So we can take the deviation of the asymptotic 
exponential from the binomial (Poissonian) distribution 
as a quantifier of the short time memory effect. Let 

g(T) = goe^ T , 

be the exponential adjusted to the asymptotic part of the 
distribution. As shown in section^ the binomial distri- 
bution results in 7 = — ln(l — fi), that reduces to 7 = fi in 
the Poissonian case. Since we are specially interested in 
small intervals, we will use the relative deviation from the 
Poisson statistics — 1 = j(T) — 1 as the quantifier 

of the memory effect. 

In FIG. |H1 it is shown the dependence of the memory 
effect with the position of the interval. We note that the 
effect occurs for virtually all recurrence intervals. Most 
of the intervals deviate from the exponential positively 
but for recurrence intervals that contain periodic orbits 
of small periods the asymptotic exponential is negatively 
deviated when compared to the binomial distribution. 




FIG. 8 Dependence of the memory effect ( r y{T) — 1) on the 
position of the recurrence interval (X c ). (a) Results for the 
10 3 intervals with measure = 10 -3 , centered in X c , that 
covers the interval x £ [0, 1] of the logistic map. (b) histogram 
of (a), (c) amplification of (a) showing the size of the fitting 
errors. The gray line represents the expected value for the 
binomial distribution. 



To explore the dependence of the memory effect on 
the size of the interval we took, by chance, an arbitrary 
position X c = 0.582 and we vary its semi-width, S. The 
results obtained numerically are illustrated in FIG.|5J For 
the logistic map, the relation between the measure, /x(/), 
and the semi-width, <5, is easily found out integrating the 
distribution (|rTi|) 



X c +8 



1 



: d,X 



x c -8 ir^/{x{l-x)) (19) 
arcsin(A c + 6) — arcsin(A c — 5) . 
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FIG. 9 Dependence of the memory effect (j{T) — 1) on the 
size of the recurrence interval (8). The position of the interval 
was chosen randomly as X c = 0.582. The dashed and solid 
lines are the expected results for the Poisson and binomial 
distribution respectively, (b) and (c) are amplifications of (a) 
for small values of 8 and show the convergence to the binomial 
statistics for 8 < 10~ 4 in the fitting precision. 

With the relation l|19f) and remembering that for the 
binomial distribution 7 = — ln(l — fJ,) we obtain the solid 
line in FIG. El From FIG. El we see that for great values 
of 8 the results deviate from both the Poissonian (dashed 
line) and binomial-like distributions. The convergence of 
the numerical results of the logistic map to the Poisso- 
nian statistics occurs in the limit of small interval, co- 
herently with what is found in the literature This 
convergence occurs clearly only for 8 < 10~ 4 . Consider- 
ing the error bars obtained by the numerical fitting, these 
are much smaller intervals than the ones for which the 
convergence of the binomial-like to Poisson takes place. 
Considering our previous discussions, the convergence of 
the numerical results to the Poissonian statistics occurs 
only when the short time memory effect becomes negli- 
gible. 



V. CONCLUSIONS 

Based on a simple Bernoulli trials problem, we obtain 
a binomial-like distribution for the r-th recurrence time 
statistics of a generic interval. This distribution depends 
only on the measure of the recurrence interval (fJ-(I)) 
that, when it is sufficiently small, turns out the statistic 
to the usual Poissonian distribution for the first Poincare 
recurrence time. 

The information related to the dynamical properties 
appears in the deviation from these distributions as, for 
example, the power-law behavior for large recurrence 
times studied in references lfT(i fl3j) . In this article, we 



discuss a kind of deviation that appears for finite size 
intervals. In this case the short recurrence times are af- 
fected by a kind of memory of the chaotic systems. We 
show that the origin of this short time memory effect is 
due to the existence of unstable periodic orbits inside of 
the finite size recurrence interval. The analytical method 
of calculating the recurrence probabilities, described in 
section IIV.AI shows how these periodic orbits changes 
the distribution of short recurrence times. Furthermore, 
our observed deviations for the specific systems studied in 
this work are in agreement with bounds predicted for gen- 
eral classes of systems in previous works (fl8t fl9t l20t l2lf) . 

The exponential growth of the number of periodic un- 
stable orbits in chaotic dynamical systems restates the 
condition of independence between recurrences for large 
times, resulting in an exponential-like behavior of the 
recurrence time distribution after the decay of mem- 
ory. Imposing the normalization condition and the Kac's 
lemma, we are able to make an analytical fitting to the 
memoryless part of the distribution. This fitting illus- 
trates how the short recurrence times, which are affected 
by the memory, modify the whole distribution. Since just 
a few points deviate from an exponential, a histogram bin 
larger than one may lead to the wrong conclusion that 
the distribution is Poissonian. What results in a con- 
tradiction: the inverse of the false Poissonian exponent 
is different from the mean recurrence time obtained from 
the recurrence series. For example, for the data of FIG. 
this difference is of 9.3%. We believe that this alert has 
to be taken into account when the RT statistics is calcu- 
lated. We would like to stress that our results indicate 
that the correct normalization of a scries of recurrence 
times should consider an extra factor, besides the mean 
recurrence time, in order to properly obtain the asymp- 
totic exponential-one law. 

We emphasize that the memory effect, discussed in this 
article, applies for any recurrence interval of a general 
chaotic dynamical system. For small intervals it may not 
be relevant because it turns smaller than the numeric 
precision. When one is calculating RT numerically, it 
must be checked that the interval is sufficiently small 
that this regime is already achieved. 
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